Coulomb oscillations in three-layer graphene 
nanostructures 



J. Giittinger, C. Stampfer, F. Molitor, D. Graf, T. Ihn, and 
K. Ensslin 

Solid State Physics Laboratory, ETH Zurich, 8093 Zurich, Switzerland 
E-mail: guettinj@ph.ys . ethz . ch 

Abstract. We present transport measurements on a tunable three-layer graphene 
single electron transistor (SET). The device consists of an etched three-layer graphene 
flake with two narrow constrictions separating the island from source and drain 
contacts. Three lateral graphene gates are used to electrostatically tune the device. An 
individual three-layer graphene constriction has been investigated separately showing a 
transport gap near the charge neutrality point. The graphene tunneling barriers show 
a strongly nonmonotonic coupling as function of gate voltage indicating the presence 
of localized states in the constrictions. We show Coulomb oscillations and Coulomb 
diamond measurements proving the functionality of the graphene SET. A charging 
energy of w 0.6 meV is extracted. 
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1. Introduction 

Carbon materials, such as carbon nanotubes and graphene have attracted increasing 
interest in the past decades, which is mainly due to their unique electronic properties. 
Strong suppression of electron backscattering in carbon nanotubes P and graphene [21 [3] 
make both materials interesting for future high mobility nanoelectronic applications [U 
[5]. Their low atomic weight and the low nuclear spin concentration, arising from the 
~ 99% natural abundance of 12 C are good premises for having weak spin orbit and 
hyperfine couplings. These make carbon nanomaterials also promising candidates for 
future spintronic devices 0, [7] and spin-qubit based quantum computation [81 191 fTUl [TT] . 
It has been shown recently that nanotubes exhibit a topologically induced spin-orbit 
coupling [12], which is directly related to their cylindrical shape [131 El- in graphene 
and few-layer graphene such flux accumulating (circumferential) trajectories should not 
be present, leading to a legitimate hope for much weaker spin-orbit interaction and 
thus possible applications for spin-based quantum information processing. However, 
graphene and few-layer graphene quantum devices are still in their infancy since it is 
hard to confine carriers in these semi-metallic materials using electrostatic potentials. 

Here we report on Coulomb oscillations and Coulomb diamond measurements on 
an etched and fully tunable three-layer graphene single electron transistor. Single 
electron transistors (SETs) consist of a small island connected via tunneling barriers 
to source and drain contacts [To] . First few-layer graphene SETs have been formed 
by Schottky barrier contacts on graphitic flakes [IE] and just very recently etched 
single-layer graphene structures [171 [18] have been fabricated to demonstrate Coulomb 
blockade. 

2. Device and fabrication 

The investigated SET device is shown in Figs. l(a)-(c) and consists of a three-layer 
graphene structure. Two ~ 60 nm wide constrictions separate the graphene island from 
the source (S) and drain (D) contacts. The two constrictions are separated by about 
1 //m, while the area of the island is 0.1 /im 2 . Two lateral side gates SGI and SG2 
allow to change the three-layer graphene barriers electrostatically and independently. 
The potential on the island can be separately changed by an additional side gate denoted 
as plunger gate (PG). The highly doped Si substrate is used as a back gate (BG) giving 
control over the overall Fermi level. Note that on the same flake also an additional single 
constriction has been fabricated (see inset in Fig. 2a). 

The devices have been fabricated by mechanical exfoliation of graphite flakes on Si 
substrates covered with 295 nm thick Si02 as described in Ref. [19]. The individual 
graphite flakes were patterned by electron-beam (e-beam) lithography using 90 nm 
PMMA as resist and a subsequent reactive ion etching step by an A1/O2 plasma (9:1). 
Fig. 1(b) shows a scanning force microscope (SFM) image of the etched three- layer 
graphene flake (bright area). A second e-beam lithography step followed by metalization 
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Figure 1. (color online) (a) Schematic illustration of the single electron transistor 
device, (b) SFM image of the device after reactive ion etching. The inset shows a 
cross section along the dashed line, (c) SFM image of the final device with metal 
contacts. The dotted lines indicate the circumference of the graphene structures, (d) 
Confocal Raman spectra recorded on the island after fabrication as highlighted in the 
Raman image shown as inset. The laser spot size is approx. 400 nm. (e) The 2D line 
can be approximated by two Lorentzians. For more details see text. 



and lift-off is used to place 2 nm Ti and 50 nm Au electrodes as shown in Fig. 1(c). 

Confocal Raman spectroscopy measurements [20] have been used to determine the 
thickness of the graphitic flake, i.e. the number of graphene layers. A Raman spectrum 
recorded at the center of the three-layer graphene island is shown in Fig. 1(d). This 
spot is marked by a circle in the Raman image [21] of the device which is shown as an 
inset in Fig. 1(d). In this Raman image, white BXCclS £1X6 attributed to the silicon oxide, 
bright (yellow) areas to the contacts and dark (blue) areas to the three-layer graphene. 
Apart from the elevated background due to nearby metal contacts the spectrum shows 
pronounced G and 2D lines typical for sp 2 graphitic materials. The defect induced D line 
arises from the edges of the flake inside the area of the laser spot, which has a diameter 
of about 400 nm. It is known from earlier experiments [201 122| 123] that the lineshape 
of the 2D peak and the intensity ratio of G/2D provides direct insight into the number 
of graphene layers of the investigated flake. Thus the 2D line is analyzed in more detail 
as shown in Fig. 1(e). According to Ref. [20] we can fit the 2D peak either by 1, 2 or 4 
Lorentzian(s), just depending on the number of graphene layers. It turns out that the 
measured spectra can be best fitted by the sum (straight line) of two Lorentzians [dashed 
lines in Fig. 1(e)]. The center of the two Lorentzians are offset by Auj = 25.0 ±0.5 cm -1 . 
In addition, the ratio between the integrated intensity of the G and the 2D line is 
~ 0.57. Both measures, including the relative height of the two Lorentzians provide 
strong support for having three-layer graphene. We can exclude bilayer graphene, since 
for two layers the integrated intensity ratio is G/2D = 0.38 ± 0.02 [20], but there is 
still a small chance of having four layer graphene. For four-layer (six-layer) graphene 
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Figure 2. (color online) High bias source-drain current measurements. (a) 
Conductance as a function of bias and back gate voltage for an individual three- 
layer graphene constriction (« 65 nm wide). A SFM image of the structure is 
shown as inset. At low bias the opening of a transport gap is observed around 
Vbg ~ —17 V. (b) and (c) show back gate characteristics for different side gate 
voltage configurations. In (b) VsG2 = and the influence of SGI is investigated 
for V SG i = 10 V (line), V (dotted), and - 10 V (dashed line). In (c) V S gi = and 
the influence of SG2 is investigated for Vsg2 = 10 V (line), V (dotted), and — 10 V 
(dashed line). A difference of the local doping of the two constrictions can be observed. 
For more information see text. 



the two Lorentzians are offset by Au = 26 ± 2 cm" 1 (Au = 28 ± 2 cm ) and the 
integrated intensity ratio is G/2D = 0.63 ± 0.08 (G/2D = 0.7 ± 0.1) [20]. It is found 
that our mesurements fit best to the value of three layers where the following values 
have been reported Auj = 25.4 ± 1.5 cm -1 and G/2D = 0.53 ± 0.05. This is also in 
good agreement with the step height of ~ 1.5 nm taken from the SFM data [inset and 
dashed line in Fig. 1(b)]. However, four- layer graphene can not be completely excluded 
here. In the following we refer to our sample to consist of three layers of graphene. The 
overall conclusions would not change if the sample had indeed four layers of graphene. 

If not stated explicitly, measurements were performed at 2 K in a variable 
temperature cryostat by applying a symmetric source-drain bias voltage Vbias (DC 
and small superimposed AC component) and measuring the source-drain current. The 
samples have been heated up in vacuum to 135° C for 12 h before cool down to eliminate 
undesired atoms on the sample surface as much as possible. 

3. Results and discussion 

We first discuss measurements on a single constriction. This structure has been 
fabricated from the same graphene flake as the SET. The investigated constriction has 
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a width of ~ 65 nm, as shown in the inset of Fig. 2(a). Differential conductance 
measurements as function of Vbias and Vbg (i- e -> Fermi level) reveal the presence of a 
transport gap around the charge neutrality point (at Vbg ~ —17 V, see white arrow). 
This is indicated by (i) conductance suppressions near the charge neutrality point, 
including strong fluctuations for low bias, which might be due to resonances in the 
three-layer graphene constriction [see dark regions in Fig. 2(a)] and (ii) the significant 
conductance non-linearity as a function of increasing bias. This result is very similar 
to what has been measured in single-layer graphene nanoribbons [2U [25], where it has 
been argued that a nanoribbon-width dependent effective energy gap (i.e. transport gap) 
dominates the transport near the charge neutrality point [26]. For the measured three- 
layer graphene constriction we observe the onset of a transport gap which is smaller than 
the thermal energy (« k-^T = 0.2 meV), since we do not observe any gap induced pinch 
off at small bias voltages. The gap is therefore significantly smaller than the ~ 4 meV, 
which has been reported for 65 nm wide single-layer graphene constrictions [26], [24]. 
However, most importantly, three-layer graphene constrictions exhibit a transport gap 
which can be used to form tunneling barriers for defining a three-layer graphene island, 
very much like in single- layer graphene [T71 [18] . 

Measurements on the SET are performed first in the high bias regime (Vbi as = 
100 mV), where transport is not suppressed by the two constrictions. This allows us 
to investigate different regimes in the back- and side gate parameter range. Figs. E](b) 
and 2(c) show the (two-point) source-drain resistance under the influence of different 
side gate potentials [Vsgi is stepped in Fig. 2(b) and Vsg2 in Fig. 2(c)]. The dotted 
traces show measurements where the side gate voltages have been set to zero. The 
dashed and the solid lines correspond to measurements where negative (Vsgi, 2 = -10 V) 
and positive (Vsgi, 2 = 10 V) side gate potentials have been applied respectively. The 
resistance shows a peak around Vbg — —15 V, which we identify as the conductance 
minimum at the overall charge neutrality point of the significantly n-doped sample. For 
increasing or decreasing back gate voltage the resistance decreases. This can be well 
explained by the (linear) carrier density increase as function of the back gate voltage 
when moving away from the charge neutrality point [19]. By applying different side 
gate potentials the peak height, width, and position change. For example, setting Vsgi 
= 10 V and sweeping Vbg the resistance peak becomes higher and narrower, while it is 
less pronounced and broader for Vsgi = -10 V [see Fig. 2(b)]. The opposite behavior 
is observed for applying ±10 V to SG2, as shown in Fig. 2(c). In both cases a positive 
side gate voltage leads to a down shift in back gate voltage, whereas a negative voltage 
leads to an up shift. 

These transport characteristics can be well explained by assuming (i) that the 
transport in this regime is dominated by the two constrictions 1 and 2 and (ii) that 
the two constrictions are differently doped, i.e., they exhibit two different (local) charge 
neutrality points. Here constriction 2 is slightly more n-doped than constriction 1. 
According to Ref. [27] graphene side gates work well for locally changing the carrier 
density, i.e., for locally shifing the charge neutrality point. Moreover, we assume 
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(will be shown below) that the crosstalk of SGI and SG2 on constrictions 2 and 1 is 
negligible. Therefore, a positive potential on SGI reduces the doping inbalance between 
constriction 1 and 2 [see Fig. 2(b)]. This increases the sample homogeneity and results 
in a high and narrow resistance peak. In contrast, by applying Vsgi = -10 V the doping 
difference increases and the resistance peak significantly shrinks and broadens. The 
measurements shown in Fig. 2(c) can be explained similarly. Additionally, we can now 
extract the relative lever ams of the side gates to the back gate with respect to the 
transport dominating constrictions given by AVsgi,2/AVbg ~2.8. From now on SGI is 
operated in a positive and SG2 in a more negative voltage regime in order to match the 
doping levels in constrictions 1 and 2, which might lead to a more symmetric coupling 
of the island to source and drain. 




Side Gate 2 Vsg2 (V) 

Figure 3. Source-drain current plotted as function of the two side gates voltages 
Vsgi and VsG2 for constant bias voltage (Vbias = 200 //V). Measurements are taken 
at V BG = -13.47 V and V PG = V. 

We now discuss the low bias transport properties. Fig. [3] shows a measurement 
of the source-drain current as function of both side gates Vsgi and Vsg2 performed 
at Vbias = 200 /iV and Vbg = —13.47 V. We observe sequences of horizontal and 
vertical stripes of suppressed current and current resonances. Their direction in the 
Vsgi-Vsg2 plane indicates that their physical origin has to be found within constriction 1 
(vertical stripes) or constriction 2 (horizontal stripes). The current exhibits even finer 
equidistantly spaced resonances which are almost equally well tuned by both side gates. 
We therefore attribute these resonances to states localized on the island between the 
barriers. It will be shown below that these resonances are Coulomb oscillations due to 
charging of the three-layer graphene island. The overall behavior is very similar to what 
has been observed in a single layer graphene SET |17j . 

Coulomb oscillations are further investigated by modulating the plunger gate 
voltage V"pg and simultaneously compensating its influence on the constrictions by the 
side gates respectively. Fig. 4 shows the current as a function of Vpc- Here SGI has 
been swept simultaneously following Vsgi = 6 V — 0.27Vpc- The plunger gate induced 
background modulation due to resonances in constriction 2 is negligible within this Vpg 
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Figure 4. (color online) Coulomb oscillations [Fig.(b)] as a function of the plunger 
gate voltage Vpg- Although the background is strongly modulated by constriction 
resonances [see Fig. (a)] an almost constant peak spacing is observed over more than 
150 oscillations [see inset in Fig. (b)]. Measurements are taken at Vsgi = 6 V— 0.27VpGj 
V SG 2 = -2.37 V, Vbg = -10.51 V and V hias = 50 aV. 



range and has therefore not been compensated. For these measurements Vbi as = 50 /xV 
has been applied and the back gate voltage has been set to Vbg — —10.51 V. Thus 
the Fermi energy in the source and drain contacts lies within the conduction band. 
In accordance with the measurement in Fig. 3 the current shows coarse and fine 
modulations. Again the larger oscillations with characteristic Vpg spacings of a few 
volts [see e.g. Fig. 4(a)] are attributed to transmission resonances in the constrictions 
while the fine current modulations at a voltage scale of around 30 mV are Coulomb 
oscillations, as shown in Fig. 4(b). The elevated background at the left side in Fig. 4(b) 
is due to resonances in the constrictions. The inset in Fig. 4(a) shows the spacing 
of the Coulomb oscillations AVp P as function of the plunger gate voltage, which has 
been swept over more than 150 periods. The gray marked region corresponds to the 
oscillations shown in Fig. 4(b). The mean Coulomb peak spacing is AV PP = 30.1 ±2 mV 
and part of the observed broadening might be due to the underlying modulation of the 
transmission through the narrow constriction [see correlation between Fig. 4(a) and 
inset in Fig. 4(b)]. The essentially constant peak spacing indicates that the three-layer 
graphene transistor compared to the single layer graphene SET [17] behaves much more 
like a metallic SET [28|. One also needs to take into account that the three-layer device 
investigated here is larger than the single layer device presented in Ref. [UJ. Therefore 
the single-level spacing which could give rise to peak spacing fluctuations is also smaller. 
Nevertheless the data in the inset of Fig. 4(b) resembles pretty much observations on 
metallic SET with the additional feature of superposed constriction resonances. Like 
for the single-layer graphene SET, some peak spacing fluctuations can be seen in Fig. 3 
indicating that weak inhomogeneities exist within the dot including its edges. 

Corresponding Coulomb diamond measurements [TH], i-e., measurements of the 
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Figure 5. (color online) Coulomb diamond measurements. From the size of the 
Coulomb diamonds we estimate a value for the charging energy of around 0.6 meV. 
The DC bias is modulated with a 50^V AC component allowing for a direct lock-in 
amplifier measurement of the differential conductance. The measurement is performed 
with Vsgi = 6 V - 0.27 • V PG , V SG2 = -2.37 V and V BG = -10.51 V. 



differential conductance (Gdiff = dI/dVbia. s ) as a function of bias voltage Vbias and 
plunger gate voltage Vpg are shown in Fig. 5. For this measurement the same gate 
voltage configuration as in Fig. 4 has been used and an AC modulation of 50 /j,V has been 
superimposed on to the DC bias. The differential conductance is plotted logarithmically 
as function of the plunger gate and bias voltage. From this measurement the charging 
energy is estimated to be Eq « 0.6 meV. This corresponds to a total capacitance of 
Cs = e 2 / 1 Eq 271 aF. The electrostatic coupling capacitances of the different gates 
to the island are Cpc ~ 7.1 aF, Csgi ~ 8.5 aF, C*sg2 ~ 7.1 aF and Cbg ~ 28.45 aF. 
The capacitance of the island to the back gate Cbg can be compared with a plate 
capacitor model leading to a capacitance C = e eA/d « 12 aF. The difference to the 
capacitance obtained from the measurement can be explained by edge effects which 
are not accounted for in the simple plate capacitor model. These measurements are in 
accordance to what has been observed for a single-layer graphene SET [17]. There the 
island area is significantly smaller, leading to a larger charging energy («3.5 meV) and 
a larger discrepancy between the back gate capacitance and the plate capacitor model 
is observed due to the even more enhanced edge effects of the smaller island. 

Finally, we investigate the temperature dependence of the Coulomb oscillations 
and the background (i.e., constriction) resonances in two different plunger gate regimes. 
Within the first regime the background is strongly elevated, whereas in the second regime 
the background is strongly suppressed since we are between two constriction resonances. 
The Coulomb oscillations on top of a constriction resonance are plotted for different 
temperatures in Fig. 6(a). Here the background has been subtracted and the traces are 
vertically offset by 10 pA for clarity. In Fig. 6(b) the data are presented in the same 
way for the "off-resonance" regime in the barrier transmission with a spacing of 3 pA 
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Figure 6. (color online) Temperature dependence of Coulomb oscillations (COs) 
in two different regimes, (a) COs sitting on top of a constriction resonance ("on- 
constriction-resonance" ) . (b) COs in a regime where the current is strongly supressed 
by the constrictions ("off-constriction-resonance"), (c) Current measured at the 
plunger gate voltages indicated by the dashed lines in (a) and (b) for both regimes, (d) 
Temperature dependence of the COs peak to peak amplitude. Measurements taken 
at V SG i = 5.067 V - 0.27(^ PG - 3.437 V), V SG 2 = -1.588 V, V^bg = -10.51 V and 
F bias = 200 iN. 

between the traces. The change of the background current is shown in Fig. 6(c), where 
the current is plotted for two fixed PG voltages indicated by the vertical dashed lines in 
Fig. 6(a) and Fig. 6(b). On top of the constriction resonance (circles) the background 
current increases linearly with temperature, which does not hold for the regime with 
suppressed transmission ("off-constriction-resonance", triangles). There, the nonlinear 
current increase is attributed to raising and subsequent broadening of constriction 
resonances lifting the overall background current. Fig. 6(d) shows the averaged 
peak to peak Coulomb oscillation amplitudes. While the "on-constriction-resonance" 
amplitudes (circles) of the Coulomb oscillations are in general decreasing with increasing 
temperature the amplitudes of the "off-constriction- resonance" oscillations (triangles) 
are mainly limited by the transmission of the constrictions and therefore increase due 
to enhanced transmission (elevated background) with increasing temperature. 

In conclusion, we have fabricated a tunable three-layer graphene single electron 
transistor based on an etched graphitic flake with lateral gates. Its functionality was 
demonstrated by observing clear and reproducible Coulomb oscillations. The tunneling 
barriers formed by three-layer graphene constrictions were investigated independently. 
From the Coulomb diamond measurements it was estimated that the charging energy of 
the three- layer graphene island is ^0.6 meV, which is compatible with its lithographic 
dimensions. The overall behavior of the investigated device is very much like that 
observed for a single-layer graphene single electron transistor. The almost constant 
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Coulomb peak spacing indicates the more metallic character of the three layer graphene 
SET. These results open the way to more detailed studies of future graphene and few- 
layer graphene quantum devices. 
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